Apex Exponents for Polymer— Probe Interactions 
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We consider self-avoiding polymers attached to the tip of an impenetrable probe. The scaling 
exponents 71 and 72, characterizing the number of configurations for the attachment of the polymer 
by one end, or at its midpoint, vary continuously with the tip's angle. These apex exponents are 
calculated analytically by e-expansion, and numerically by simulations in three dimensions. We find 
that when the polymer can move through the attachment point, it typically slides to one end; the 
apex exponents quantify the entropic barrier to threading the eye of the probe. 
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There has been remarkable progress in recent years 
in nanoprobing and single-molecule techniques. These 
developments have had a direct impact on biopolymer 
research producing a wealth of beautiful results on DNA 
dynamics hi, molecular motors 0, and protcin/RNA 
folding 0, Today it is possible to measure statis- 
tical properties of a single macromolecule rather than 
deducing them from experiments with solutions of many 
polymers. This naturally leads to questions regarding the 
theoretical limitations of these techniques, such as the ef- 
fects of microscopic probes on the measured properties of 
the polymer. Consider, for instance, a polymer attached 
to the apex of a cone-shaped probe (e.g. a micropipette 
or the tip of an atomic force microscope HE1). What 
is the configurational entropy for this system? Suppose 
that this probe is a microscopic needle with a hole at the 
end. How hard is it to thread a polymer through the 
needle's eye? 

Quite generally, the number of configurations 9t of a 
polymer of length A or, equivalently, of an A-step self- 
avoiding walk (SAW), behaves as |j| 

m = const x z^A 7 " 1 . (1) 

The "effective coordination number" z, depends on mi- 
croscopic details, while the exponent 7 is 'universal.' Ac- 
tually, 7 does depend on geometric constraints which in- 
fluence the polymer at all length scales. In particular, 
there are a number of results demonstrating the varia- 
tions of 7 for polymers confined by wedges in two and 
three dimensions [1 H EH EH : A SAW anchored at the 
origin and confined to a solid wedge (in 3D) or a pla- 
nar wedge (in 2D) has an angle-dependent exponent 7 
that diverges as the wedge angle vanishes. A limiting 
case which has been extensively studied, both analyt- 
ically [H E3 and numerically [H E3, is a SAW an- 
chored to an impenetrable surface, for which 7 = 7,, = 

0.70 ±0.02 [n[. 

To model the polymer-probe system, we consider a 
SAW attached to the apex (tip) of an impenetrable obsta- 




FIG. 1: Configurations of a polymer near an obstacle: (a) at- 
tached to the apex of a planar sector of angle a; (b) threaded 
through the eye of a cone with apex semi-angle /3. 

cle (needle) . To avoid introduction of an external length 
scale, we focus on obstacles of scale-invariant shape, such 
as a planar slice (sector) of angle a (Fig. ^) , or a con- 
ical needle of apex semi-angle f3 (Fig. ^p) . While both 
geometries are natural extensions of the 2D wedge, they 
are clearly different in three dimensions (and also distinct 
from the 3D wedge, which consists of two planes inter- 
secting at a line). The former excludes the polymer from 
the volume of a cone, while the latter prevents it from 
crossing the surface of a slice. Nonetheless, the result- 
ing phenomenology is rather similar. Indeed, one of the 
technical innovations of this paper is the demonstration 
that many such geometries can be treated in the same 
manner by an e — 4 — d expansion focusing on the in- 
teraction with a 2D surface. The e-expansion, as well 
as numerical simulations in 3D, shows that the exponent 
7 = 71 varies continuously with the apex opening angles 
in Fig.^ Continuously varying exponents are rather un- 
common in critical phenomena. In the present case they 
arise from the interaction of two self-similar entities, the 
polymer and the probe. 

Another variant of this problem occurs when a poly- 
mer is attached to the apex at its midpoint. This case is 
described by Eq. 0) with exponent 7 = 72. More gener- 
ally, let us denote by 9^2 (A, Ai), the number of accessible 
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configurations for a polymer attached to the apex at an 
arbitrary monomer, dividing it in two segments of lengths 
Ni and N 2 = N — Aq . If we allow the two segments to ex- 
change monomers with each other (which can be done by 
replacing a rigid attachment with a slip-ring as depicted 
in the Fig. (IJd) , then the equilibrium configurations will 
be distributed with a weight proportional to 0^2 (N, Aq). 
A natural interpolation formula as a function of a (or (3) , 
supported by the e-expansion at first order, is 



Vt 2 (N,Nx) cx AT c ^[iVi(iV-iVi)] Cl H 



(2) 



To get a feeling for this scaling relation, let us look at 
some limits: When the probe is absent, we recover Eq. (fTJl 
and c(0) = To — 1, where 70 — 1.158 describes the geomet- 
rically unconstrained SAW. If the obstacle is present but 
the two segments do not interact with each other, then 
c = and C\ = 71 — 1 . By fitting to the limits of Aq — > 
and Aq ~ N%, we find C\ = 72 — 71 and c = 27! — 72 — 1. 
Below, we estimate the exponents in Eq.(j2J) both analyt- 
ically and numerically. For now, assuming Eq. © holds, 
we see that if c\ < 0, the maximum number of configu- 
rations is realized when either Aq or N 2 equals N. This 
brings us to one of our main findings: No matter how 
small the apex angle, we find c\ < 0, i.e. the most likely 
states have Aq ~ N or N 2 — N, with an entropy barrier 
separating the two. Threading a needle is hard! 

To treat the problem analytically, we start with the 
Edwards 0] model of a self- avoiding polymer, and add 
an interaction with the obstacle. In this formulation, 
configurations of the polymer are described by r(r) £ 5R d , 
where r measures the position along the chain, and are 
weighted according to the energy |24j 

U = - I r 2 dr + -^- dr dT'6[r( T ) - r(r')] 
2 Jo 2 Jo J Q 

+ go [ d 2 R [ rfr5[r(r) - R]. (3) 

J M JO 

The self-avoiding interaction is replaced by a 'soft' repul- 
sion of strength Vq. In the same spirit, the impenetrable 
obstacle is replaced with a soft repulsion of magnitude 
go- The key observation is that in 3D the polymer can 
only sense the exterior of an impenetrable obstacle, and 
will not care if its interior is hollow. In generalizing to 
e?-dimensions, we keep the dimensions of the now soft- 
ened exterior manifold (indicated by R £ M) as two. 
The advantage of this choice is that both 170 and vq 
have the same bare dimensions, and in a perturbative 
scheme simultaneously become relevant in d < 4. We 
then analyze the model using a renormalization group 
(RG) scheme [HEH which is a modification of the con- 
formation space RG GI0]|. The scaling exponents are 
calculated using dimensional regularization in d = 4 — e 
dimensions to order O(e). 

It is customary to define non-dimensionalized (bare) 
coupling constants i>o — voL e , go = goL e at a length scale 




FIG. 2: Diagrams contributing to renormalization of g to 
second order (a-c); to Z in first order (d,e) at the apex of a 
slice; and to -Z2 in first order (f) at the eye of a conic needle. 



L. We also define the renormalized coupling constants 
v = Z~ l vo and g = Z~ go, where the renormalization 
constants Z v , Z g are calculated perturbatively as series in 
v and g. Diagrams contributing to the renormalization of 
g are shown in Figs.|2K - c, and involve both interactions 
of the polymer with the obstacle and with itself. The 
leading singularity in 1 /e comes from short distances and 
therefore the leading correction to g does not depend on 
the overall shape of the obstacle. This is true as long as it 
is possible to draw a finite circle around points away from 
the apex (see Fig. |2t) , and is also why the interaction 
with the obstacle becomes irrelevant in the degenerate 
limits of a — (3 — 0. The self-interaction is uninfluenced 
by the obstacle, and the renormalization constant Z v is 
the same as in the unattached polymer |l6j . To first 
order in e, the nontrivial fixed point of the RG flow is 
found to be (v*,g*) = (7r 2 e/2, 37re/4). The fixed point 
thus depends only on the dimension of the constraining 
manifold, but not on its shape [ItI H^ |. However, the 
number of accessible configurations does depend on the 
exact geometry as described below. 

Consider first-order corrections to the partition func- 
tion Z coming from the self-interaction (Fig.|2i) and the 
interaction with the slice (Fig. [5^). Combining them and 
adding relevant counterterms to eliminate poles in e, we 
obtain 



Z = 1+ 4>* 



a 5 *)ln(^)+0( e 2 ). (4) 



Comparing this with i\T7i(«)-i = 1+ (ct) — 1) In iVH , 

and substituting the fixed point values (v*,g*), we find 



71(a) = 1 + 



1 



3 a 
2^ 



+ 0(e 2 ). 



(5) 



The above treatment is easily generalized to a polymer 
attached by its midpoint. For the number of configura- 
tions, we observe that the contribution from the inter- 
action with the obstacle is doubled. Interaction between 
the two halves of the polymer, however, makes no sepa- 
rate correction and is already included. (Note that if we 
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ignore the obstacle and consider self-interactions only, we 
get a "degenerate" star polymer with two branches that is 
equivalent to a linear polymer.) Thus, for the calculation 
of 72 at order of e, we can add the separate contributions 
from self-avoidance and avoidance of the obstacle; cross- 
terms can only occur at higher orders. This enables us 
to identify the scaling exponent 

72(a)-l + ^(l-^)+0(£ 2 ). (6) 

Repeating this argument for the slip-ring geometry, we 
find 

m 2 oc N e / S [Ni (N - Ni)]- 3 ^^, (7) 

which confirms the Ansatz in Eq. to first order in e. 
Note that in 3D we must have 72(271) = 2j s — 1. This 
equality does not hold in the e-expansion. The reason 
is that in 3D, a complete plane prevents two polymers 
on its opposite sides from interacting with each other, 
whereas in 4D it does not. 

It is straightforward to extend the above formalism to 
obstacles of different shapes, such as the conical manifold 
with apex angle (3 (Fig. ^3). In counting the number of 
configurations we obtain a result similar to Eq. Q , with 
a replaced by 2irs'm(3. Since the fixed point location is 
the same as before, this substitution in Eqs. Il- r >l(j[| gives 

7l conc (/3) = l + ^(l-3sin/3) + 0(e 2 ), (8) 

Wf nc (/?) oc 7V e / 8 [A^ 1 (^-7V 1 )]-( 3 / 8 ) esin ' 3 . (9) 

The difference between the two geometries is thus merely 
quantitative. As another point of caution, we note that 
within our formalism, the polymer is free to occupy either 
side of the hollow cone- the partition sum is dominated 
by the arrangement with the largest number of configura- 
tions. Thus the result for 71 (/3) is valid only for (3 < tt/2, 
with a possibly more severe restriction for 72. 

The earlier discussion of 'threading a needle' illus- 
trates the essence of the method of "entropic competi- 
tion" [H |3 which we employ to numerically estimate 
the exponents 74(a) in 3D. The underlying idea is to sam- 
ple the ensemble of different configurations of two poly- 
mer segments which can exchange monomers and thus 
"compete entropically" (e.g. in a slip-ring geometry). 
To calculate 71(a), we prevent the two segments from 
interacting with each other. The number of configura- 
tions for given N and N2 = N — Ni is then 

3ti oc [Nx(N - Ni)]^- 1 , (10) 

so that the resulting histogram for N\ allows us to calcu- 
late the exponent 71(a). In our simulations we consider 
sectors with angles a — kir/8, where k = 0, 1, 2, . . . , 16. 
We work on a cubic lattice, with the excluded sector in 
the x — y plane. Possible Monte Carlo (MC) moves in- 
clude attempts to remove one monomer from the free end 
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FIG. 3: The probability distributions p(Ni) for two non- 
interacting segments of lengths Ni and N — Ni attached to 
the apex of a planar slice for diffrent values of angle a. The 
curves are the result of 10 9 MC steps for N = 2000. 

of a randomly chosen polymer segment and adding it to 
the free end of the other segment; both segments also 
undergo random configuration changes via pivoting |2l| . 
Figure illustrates the dramatic effect of the angle a on 
p(Ni), the probability distribution function (PDF) for 
the segment length N\. For small a, the distribution is 
peaked at the center while for a bigger than a critical 
value a c , the maximum of the PDF moves to the sides. 
The numerical data from entropic competition suggest 
a c ~ 57r/8, which is not too far from the first order e— 
expansion result of a c — 2ir/3 in Eq. J5|. 

For the purpose of calculating 72(a), it is necessary 
to include interactions between the segments. Figure 0] 
depicts variations of the exponents 71,2(0) obtained by 
fitting histograms from entropic competition, such as in 
Fig. 01 to power-laws as in Eqs. I|10|) and J5J. 

It is instructive to compare the results of entropic 
competition with those of a more established proce- 
dure, such as dimerization [2^, Ej|. The latter is a 
quite efficient method [2l|, in which an A-step SAW 
is created by generating two (A/2)-step SAWs and at- 
tempting to concatenate them. We generated SAWs for 
N = 16, 32, • • • , 2048, and by attempting to attach them 
to the end point of an appropriate sector, measured a suc- 
cess probability p^. Let us indicate the number of SAWs 
not attached to the sector by AqZ N N 70 ^ 1 , and those at- 
tached to the sector either (1) by their ends, or (2) by 
their mid-point as AiZ N N 1 ^ ^ 1 (i = 1, 2 corresponds to 
the notation introduced earlier). Then, the ratio between 
the number of configurations, p^ = (Ai/A^N 1 ^ 01 ^ 1 ", 
represents the probability to attach an A^-step polymer 
to a sector with angle a. Fitting a power law to this 
ratio thus provides a means of estimating the exponent 
difference A7; =70 — 74. 

Using the dimerization method we generated M = 10 6 
SAWs; for each we checked whether there is an intersec- 
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FIG. 4: Extrapolated values of the exponents A71 = 70 — 71 
(circles) and A72 =70 — 72 (diamonds) as a function of sec- 
tor angle a from 'entropic competition' (open symbols), and 
dimerization (full symbols). Error bars represent statistical 
uncertainties of individual estimates of the exponents, as well 
as the uncertainty in the extrapolation N — > 00. 

tion with sectors at all tested values of a [2^] . The error 
in the measurement of p^ is \/pn(1 — Pn)/M ■ Since p n 
decreases with increasing angle of the sector, the relative 
error in p^, and an error in the estimates of A7^, in- 
creases with increasing a. Nevertheless, we were able to 
obtain reasonable estimates of the exponent for all values 
of a, as shown in Fig. 0] (full symbols). 

The two numerical approaches are in very good agree- 
ment; error bars for 'entropic competition' results being 
even smaller. For a = 0, our results deviate from zero 
beyond the statistical error range. We believe this devi- 
ation to be a finite size effect, due to discreteness of the 
lattice. As a check, we estimated A7i : 2 when the obsta- 
cle consists of the positive £-axis. While asymptotically 
such a situation corresponds to a = 0, and should lead to 
A7 l = 0, we obtained A71 = 0.02 and A72 = 0.05. Fol- 
ey = 2tt, we expect to have A71 = 70 — 7 S « 0.46, and 
A72 = 70 — + 1 « 0.76; our results are quite close to 
these estimates. 

In summary, we consider configurations of a polymer 
attached to the apex of a self-similar probe (at least on 
the scale of the polymer size). The geometric constraints 
imposed by the impenetrable probe lead to exponents 
continuously with the apex angle. Two such exponents 
are associated with attachment of the polymer by one 
end or by a mid-point. Together, they determine if a 
mobile attachment point is likely to be in the middle or 
slide to one side. These apex exponents are obtained 
analytically by an e = 4 — d expansion, and through in- 
dependent numerical schemes in d — 3. The e-expansion 
takes advantage of the marginality of interactions of a 



polymer with a two dimensional manifold in four dimen- 
sions, and can be applied to a variety of shapes. The 
numerical method of 'entropic competition' is shown to 
be a powerful tool in this context, comparable to or bet- 
ter than the more standard dimerization approach. The 
numerical and analytical results are in agreement, and 
indicate the presence of an entropic barrier that favors 
attachment of the polymer to the apex at its end. It 
would be interesting to see if these predictions can be 
probed by single molecule experiments. 
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